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Abstract 

This article is a survey concerning the state-of-the-art mathematical 
theory of the Euler equations of incompressible homogenous ideal fluid. 
Emphasis is put on the different types of emerging instability, and how 
they may be related to the description of turbulence. 



1 Introduction 

This contribution is mostly devoted to the time dependent analysis of the 2d 
and 3d Euler equations 

d t u + V • (u ® u) + Vp = , V • u = , (1) 

of incompressible homogenous ideal fluid. We intend to connect several known 
(and maybe less known) points of view concerning this very classical problem. 
Furthermore, we will investigate the conditions under which one can consider 
the above problem as the limit of the incompressible Navier-Stokes equations: 

d t u v + V • (u u ® u v ) - v/\u y + Vp y = , V • u v = , (2) 

when the viscosity v — * 0, i.e., as the Reynolds number goes to infinity. 

At the macroscopic level the Reynolds number, Re, corresponds to the ratio 
of the strength of the nonlinear effects and the strength of the linear viscous 
effects. Therefore, with the introduction of a characteristic velocity, U, and a 
characteristic length scale, L, of the flow one has the dimensionless parameter: 

Re = — . (3) 
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With the introduction of the characteristic time scale T = L/U and the dimen- 
sionless variables: f 

x = — , t = — and u = — , 
the Navier-Stokes equations $2$ take the non-dimensional form: 

d t v! + V x > ■ (u' ® u') - -p-AaV + V x >p' = , V • u' = . (4) 
Re 

These are the equations to be considered in the sequel, omitting the ' and 
returning to the notation v for Re~ x . 

In the presence of physical boundary the problems (fTJ) and ([2|) will be con- 
sidered in the open domain ft C R d ,d=2, d = 3 , with a piecewise smooth 
boundary dil . 

There are several good reason to focus at present on the "mathematical 
analysis" of the Euler equations rather than on the Navier-Stokes equations. 

1. Turbulence applications involving the Navier-Stokes equations ^ often 
correspond to very large Reynolds numbers; and a theorem which is valid for 
any finite, but very large, Reynolds number is expected to be compatible with 
results concerning infinite Reynolds numbers. In fact, this is the case when 
Re = oo which drives other results and we will give several examples of this 
fact. 

2. Many nontrivial and sharp results obtained for the incompressible Navier- 
Stokes equations rely on the smoothing effect of the Laplacian, with viscosity 
v > 0, and on the invariance of the set of solutions under the scaling: 

u(x,t) i-> Aw(Ax, X 2 t) . (5) 

However, simple examples with the same scalings, but without a conservation 
law of energy may exhibit very different behavior concerning regularity and 
stability. 

1. With <j) being a scalar function, the viscous Hamilton- Jacobi type or 
Burgers equation 

d t <t>- fA0+ i|V0| 2 = Oin Q x R+ , (6) 

4>(x,t) = for x€dQ, and <f>{-,0) = (f> (-) € L°°(fJ) , 

has (because of the maximum principle) a global smooth solution, for v > 0. 
However, for v = 0, it is well known that certain solutions of the inviscid Burgers 
equation ([6]) will become singular (with shocks) in finite time. 

2. Denote by |V| the square root of the operator —A, defined in £1 with 
Dirichlet homogeneous boundary conditions. Consider the solution u(x,t) of 
the equation 

d t u- vAu+ ^|V|(ii 2 ) = in Q x R+ , (7) 
u(x,t) = for x e dn, and w(-,0) = w o (0 € L°°(Q) . (8) 
Then one has the following proposition. 
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Proposition 1.1 Assume that the initial data uq satisfies the relation: 

u {x)(j)i(x)dx = -M < , (9) 



where 4>\{x) > denotes the first eigenfunction of the operator —A (with Dirich- 
let boundary condition), — A0i = \\4>i . Then if M is large enough, the corre- 
sponding solution u(x,t) of the system ^ J5J) blows up in a finite time. 

Proof. The L 2 scalar product of the equation ([7|) with <fri(x) gives 
— f u(x, t)<t>i(x)dx + uXi f u(x,t)(f>i(x)dx 



dt 



u(x, t) 2 cf)i(x)dx . 
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Since <f>\ (x) > then the Cauchy-Schwarz inequality implies 



u(x, t)(j>i(x)dx J — u { x -ii) 4>i{x)dxj \J (j>i(x)dx 

As a result of the above the quantity m(t) = — J n u(x, t)4>\(x)dx satisfies the 

relation: 

dm , vAi f <t>i(x)dx 9 

+ Aim > Jn — m 2 , withm(O) = M ; 

dt 2 w » 

and the conclusion of the proposition follows. 

Remark 1.1 The above example has been introduced with Q = R 3 by Montgomery- 
Smith [55] under the name "cheap Navier-Stokes equations" with the purpose 
of underlying the role of the conservation of energy (which is not present in 
the above examples ) in the Navier-Stokes dynamics. His proof shows that the 
same blow up property may appear in any space dimension for the solution of 
the "cheap hyper-viscosity equations" 

d t u + v{-A) m u + i|V|(u 2 ) = . 

On the other hand, one should observe that the above argument does not apply 
to the Kuramoto-Sivashinsky-like equations 

9 t + ^(-A)"V + aA0+i|V0| 2 = 0, (10) 

for m > 2. Without a maximum principle or without the control of some sort of 
energy the question of global existence of smooth solution, or finite time blow up 
of some solution, to the above equation is an open problem in R™, forn > 2 and 
for m > 2. However, if in ilO\) the term |V</>| 2 is replaced by |V</>| 2+7 ,7 > 
one can prove the blow up of some solutions (cf. T7j and references therein). 
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In conclusion, the above examples indicate that the conservation of some 
sort of energy, which is guaranteed by the structure of the equation, is essential 
in the analysis of the dynamics of the underlying problem. In particular, this 
very basic fact plays an essential role in the dynamics of the Euler equations. 

Taking into account the above simple examples, the rest of the paper is or- 
ganized as follows. In section 2 classical existence and regularity results for the 
time dependent Euler equations are presented. Section 3 provides more exam- 
ples concerning the pathological behavior of solutions of the Euler equations. 
The fact that the solutions of the Euler equations may exhibit oscillatory behav- 
ior implies similar behavior for the solutions of the Navier-Stokes equations, as 
the viscosity tends to zero. The existence of (or lack thereof) strong convergence 
is analyzed in section 4 with the introduction of the Reynolds stresses tensor, 
and the notion of dissipative solution. A standard and very important problem, 
for both theoretical study and applications, is the vanishing viscosity limit of so- 
lutions of the Navier-Stokes equations subject to the no-slip Dirichlct boundary 
condition, in domains with physical boundaries. Very few mathematical results 
are available for this very unstable situation. One of the most striking results 
is a theorem of Kato [38], which is presented in section 5. Section 6 is again 
devoted to the Reynolds stresses tensor. We show that with the introduction of 
the Wigner measure the notion of Reynolds stresses tensor, deduced from the 
defect in strong convergence as the viscosity tends to zero, plays the same role 
as the one originally introduced in the statistical theory of turbulence. When 
the zero viscosity limit of solutions of the Navier-Stokes equations agrees with 
the solution of the Euler equations the main difference is confined in a boundary 
layer which is described by the Prandtl equations. These equations are briefly 
described in section 7. It is also recalled how the mathematical results are in 
agreement with the instability of the physical problem. The Kelvin-Helmholtz 
problem exhibits also some basic similar instabilities, but it is in some sense 
simpler. This is explained at the end of section 7, where it is also shown that 
some recent results of [H] , [7T] and [72], on the regularity of the vortex sheet 
(interface), do contribute to the understanding of the instabilities of the original 
problem. 



2 Classical existence and regularity results 
2.1 Introduction 

The Euler equations correspond, formally, to the limit case when the viscosity 
is 0, or the Reynolds number is infinite: 

d t u + V • (u <g> u) + Vp = , V • u = , in fl . (11) 

In the presence of physical boundaries, the above system is supplemented with 
the standard, no-normal flow, boundary condition: 

u-n = on ffl, (12) 
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where n denotes the outwards normal vector to the boundary dQ. It turns out 
that the vorticity, oj — V A u, is "the basic quantity" , from both the physical 
and mathematical analysis points of view. Therefore, equations (jTTJ) and (|12|) , 
written in terms of the vorticity, are equivalent to the system: 

dtLo + u ■ Vlo = to ■ Vit in O (13) 
V-u = 0,VAu = £JinO, and u ■ n = on dfi . (14) 

That is, system (|14p fully determines u in terms of u>, which makes the above 
system "closed" . More precisely, the operator K : to i— > u defined by relation Q14p 
is a linear continuous map from C a (fl) with values in C a+1 (f2) (with a > 0); 
and from H s (£l) with values in H s+1 (il) . 

Furthermore, for 2d flows, the vorticity is perpendicular to the plane of 
motion and therefore equation (|13|) is reduced (this can also be checked directly) 
to the advection equation 

dtto + u-Vuj = 0. (15) 

The structure of the quadratic nonlinearity in (|13p has the following conse- 
quences, which are described below. We will be presenting only the essence of 
the essential arguments and not the full details of the proofs (see, e.g., [50 or 
[51] for the details). 



2.2 General results in 3d 

The short time existence of a smooth solution for the 3d incompressible Eulcr 
equations has been obtained already a long time ago, provided the initial data 
are smooth enough. To the best of our knowledge the original proof goes back 
to Lichtenstein [45] . The proof is based on a nonlinear Gronwall estimate of the 
following type: 

y> < Cyi => y(t) < ^ . (16) 

W " (l-2«7y3(0)) 2 

Therefore, the value of y(t), which represents an adequate norm of the solution, 
is finite for a finite interval of time; which depends on the size of the initial value 
of y(0), i.e. the initial data of the solution of Euler. These initial data have to 
be chosen from an appropriate space of regular enough functions. In particular, 
if we consider the solution in the Sobolev space H s , with s > | , then by taking 
the scalar product, in the Sobolev space H s , of the Euler equations with the 
solution u, and by using the appropriate Sobolev estimates we obtain: 



ld\\u\\ 2 H 



2 dt 



-(V- {u®u),u) H > <C,H^'|V«||l» <C\\uf H s. (17) 



As a result of (fTB]) and (|17[) we obtain the local, in time, existence of a smooth 
solution. 

As standard in many nonlinear time dependent problems local regularity of 
smooth (strong) solutions implies local uniqueness and local stability (i.e, con- 
tinuous dependence on initial data). Furthermore, one may exhibit a threshold 
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for this existence, uniqueness, and propagation of the regularity of the intial 
data (including analyticity Bardos and Benachour [3])- More precisely one uses 
the following theorem. 

Theorem 2.1 Beale Kato Majda fSjj Letu(t) be a solution of the 3d incom- 
pressible Euler equations which is regular for < t < T ; that is, 

for all t S [0, T] , u(t) € H s {n) , for some s > - . 

3 

Assume that 

T 

||V Au(.,t)\\ L ^dt < oo, (18) 

then u{t) can be uniquely extended up to a time T + S (5 > 0) as a smooth 
solution of the Euler equations. 

The main interest of this statement is the fact that it shows that if one starts 
with smooth initial data, then instabilities appears only if the size of the vorticity 
becomes arbitrary large. 

Remark 2.1 The Beale-Kato-Majda theorem was first proven in the whole 
space in |3]/. Extension to a periodic "box" is easy. For a bounded domain with 
the boundary condition u-n — it was established by Ferrari [26]. By combining 
arguments form [3] and [261 one can show, as in the Beale-Kato-Majda theorem, 
that the solution of 3d Euler equations, with real analytic initial data, remains 
real analytic as long as H8\) holds. 

The Beale-Kato-Majda result has been slightly improved by Kozono [40] who 
proved that on the left-hand side of (|18p. the || • l^oo-norm can be replaced by the 
norm in the BMO space. This generalization is interesting because it adapts 
harmonic analysis (or Fourier modes decomposition) techniques which is an 
important tool for the study of "turbulent" solutions; indeed, the space BMO, as 
the dual space of the Hardy space TL 1 , is well defined in the frequency (Fourier) 
space. In fact, cf. [52], BMO is the smallest space containing L°° , which is 
also invariant under the action of a zero order pseudodifferential operators. The 
idea behind the Beale-Kato-Majda theorem, and its improvement, is the fact 
the solution u of the elliptic equations (TH]) satisfies the following estimate, for 
1 < p < oo , 

UVulw-*.* < C SiP {\\u\\ W s, P + \ulw'.p) ■ (19) 

This relation could also be phrased in the context of Holder spaces C k,a , a > 0. 
We stress, however, that the estimate flT9l ceases to be true for p = oo (or 
a = 0.) This is due to the nature of the singularity, which is of the form 
\x — y\ 2 ~ d , in the kernel of the operator K and which leads (for s > d/2 + 1) to 
the estimate: 

|Vu|U- < C (Mi- log(l + \\u\\ 2 HS )) , (20) 
or sharper ||Vu|l« < C (\uj\bmo log(l + \u\%.)) . (21) 
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With z = 1 + II ufjjB and thanks to (|2ip the inequality (|17p becomes: 

^-z < C\\uj\\ bmo zlogz . 
at 

This yields: 

(i + IKt)ll^)<(i + IKo)ll^) 



2 \ / /1 il /nMl2 \e c So ll"W»BMO^ 



which proves the statement. The uniqueness of solutions can be proven along 
the same lines, as long as 

\\v(s)\\bmo 

Jo 

remains finite. 

Remark 2.2 The vorticity uj can be represented by the anti- symmetric part of 
the deformation tensor Vit . However, in the estimates \20\) or 121]) this anti- 
symmetric part, i.e. lu, can be replaced by the symmetric part of the deformation 
tensor 

S(u) = -(Vu+(Vu) t ). (22) 

Therefore, the theorems of Beale~Kato~Majda and Kozono can be rephrased in 
terms of this symmetric tensor \^0j - 

In fact, the above deformation tensor S(u) (or S(u>) when expressed in term of 
the vorticity) , plays an important role in a complementary result of Constantin, 
Fcffcrman and Majda [16] . which shows that it is mostly the variations in the 
direction of the vorticity that may produce singularities. 

Proposition 2.1 \16$ Let u, which is defined in Q = £1 x (0,T), be a smooth 
solution of the Euler equations. Introduce the quantities k\(t) and k<z{t) (which 
are well defined for t < T): 

ki(t) = sup \u(x, t) \ , 
which measures the size of the velocity, and 

fc 2 (i)=47r sup 1 _ , 

x,y(E.Q , x^£y |*** U\ 

which measures the Lipschitz regularity of the direction 

* X ' t} \u,(x,t)\ 
of the vorticity. Then under the assumptions 

(fci(t) + k 2 (t))dt < oo and / k 1 (t)k 2 (t)dt < oo , (23) 

Jo 

the solution u exists, and is as smooth as the initial data up to a time T + 5 for 
some 8 > . 
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Proof. As before we only present here the basic ideas, and for simplicity we 
will focus on the case when f2 = R 3 . First, since 

S(u) = l - {Vu + (Vu)«) (x, t) = S(u) (x,t), (24) 

we have 

i(<9 t M 2 +u- VM 2 ) = {u-Vu,u>) = {S(u))w,u>), (25) 



which gives 



!^<sup(|S( W )|)Moo. (26) 



Next, we consider only the singular part of the operator u> i— » S'(w) . The Biot- 
Savart law reproduces the velocity field from the vorticity according to the 
formula: 

, x 1 /" — v)Au)(y) , 

= 4^1 \x-y? (2?) 

For the essential part of this kernel, we introduce two smooth nonnegative radial 
functions (3( and (3^ with 

/3f + $ = 1 , p\ = for |a;| > 2<5 and /3 2 = for \x\ < 5 . (28) 

Then we have 

\S{u>)\ < I /( |^| ' £(*))(Det(^, + yU{x))(3}{\y\)\u{x + y)\ ^| + 
(^ • ^( a; ))(Det(||, £(x + j,), £(z))/3 2 (|y|))|u;0r + y)| ^| . (29) 
For the first term we use the bound 

(Det(^,e(x + 2/U(x))/?](M)) < ^\y\. (30) 

to obtain: 

(|| • £(x))(Det(^ , £(x + y), ^^(lyDk^ + i/)lpH 
< ka(t)5Moo (31) 
Next, we write the second term as 

J ( • C(x))(Det(^ , e(x + y), £(aO)$(fol))(£(s + y) ' (V, A «(* + y))) ^ 



and integrate by parts with respect to y. With the Lipschitz regularity of £ one 
has 

|V„ f • £(*))(Det( j^y + !/),£(*)) J | < Cfc 2 (i) . 
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Therefore, one has (observing that the terms coming from large values of \y\ 
and the terms coming from the derivatives of /3f (|y|) §^ ve more regular contri- 
butions.) 

| J{1L . # X ))^*-Z( X + y),Z( X ))ffi(\y\))\u(x + V)lpfl 



< 



/ |v y (|| ■m)(^t(^^ + y),m)M(\y\))^\W\\c 



<Ck2(*)|l°g(*)IMoo 
<Cfci(t)A a (*)|IogJ|. 

Finally, inserting the above estimates in (|26|) one obtains for \\u>\\oo > 1 and 
§ = \\"U 

^ Ck 2 (t)(l + fciWJIwIoo log [wloo , 

and the conclusion follows as in the case of the Beale-Kato-Majda Theorem. 

The reader is referred, for instance, to the book of Majda and Bertozzi [51)] 
and the recent review of Constantin 1141 for addition relevant material. 



2.3 About the two-dimensional case 

In 2d case the vorticity lu — \7 A u obeys the equation 

d t (VAu) + {u- V)(VAit) = 0. (32) 

This evolution equation guarantees the persistence of any LP norm (1 < p < oo) 
of the vorticity. Taking advantage of this observation Youdovitch proved in his 
remarkable paper [74j the existence, uniqueness, and global regularity for all 
solutions with initial vorticity in L°° . If the vorticity is in L v ', for 1 < p < oo, 
then one can prove the existence of weak solutions. The same results hold also 
for p = 1 and for vorticity being a finite measure with "simple" changes of sign. 
The proof is more delicate in this limit case, cf. Delort [21] and the section I7T21 
below. 



3 Pathological behavior of solutions 

Continuing with the comments of the previous section one should recall the 
following facts. 

• First, in the three-dimensional case. 

i) There is no result concerning the global, in time existence of smooth 
solution. More precisely, it is not known whether the solution of the Euler 
dynamics defined with initial velocity, say in H s , for s > | + 1, on a finite time 
interval can be extended as a regular, or even as a weak, solution for all positive 
time. 

ii) There is no result concerning the existence, even for a small time, of a 
weak solution for initial data less regular than in the above case. 



9 



iii) Due to the scaling property of the Euler equations in R 3 , the problem 
of global, in time existence, for small initial data, is equivalent to the global 
existence for all initial data and for all t G R . 

• Second, both in the 2d and the 3<i cases, the fact that a function u G 
L 2 ([0, T]; L 2 (R d )) is a weak solution, i.e., that it satisfies the following rela- 
tions in the sense of distributions 



is not enough to define it uniquely in terms of the initial data (except in 2d with 
the additional regularity assumption that V A u G L°° .) More precisely, one 
can construct, following Scheffer [64 and Shnirclman 65j, both in 2d and 3d, 
nontrivial solutions u G L 2 (R t ; L 2 (R d )) of (j3"3"j) that are of compact support in 
space and time. 

The following examples may contribute to the understanding of the underly- 
ing difficulties. First, one can exhibit (cf. Constantin [9], Gibbon and Ohkitani 
|31j . and references therein) blow up for smooth solutions, with infinite energy, 
of the 3d Euler equations. Such solutions can be constructed as follows. The 
solution u is (xi, X2) periodic on a lattice (R/LZ) 2 and is defined for all X3 G R 
according to the formula 

U = (ui(xi,X 2 ,t),U 2 (xi,X2,t),X3j(xi,X2,t)) = (u, X37) , 

which is determined by the following equations: 

To maintain the divergence free condition, it is required that 

V ■ u + -f = , 

and to enforce the Euler dynamics, it is required that 

d t (V A u) + (it • V)(V A u) = 7U 
d n + (u ■ V) 7 = - 7 2 + I(t) , 

and finally to enforce the (xi,x 2 ) periodicity it is required that 



L " J[0,L] 2 

Therefore, the scalar function 7 satisfies an integrodifferential Ricatti equation 
of the following form 



from which the proof of the blow up, including explicit nature of this blow up, 
follows. 

The above example can be considered as non-physical because the initial 
energy 



d t u + V • (u «) u) + Vp = , V • u = 0, u(x, 0) = u (x) , 



(33) 
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is infinite. On the other hand, it is instructive because it shows that the conser- 
vation of energy, in the Euler equations, may play a crucial role in the absence of 
singularity. Furthermore, an approximation of the above solution, by a family of 
finite energy solutions, would probably be possible but to the best of our knowl- 
edge this has not yet been done. Such an approximation procedure would lead 
to the idea that no uniform bound can be obtained for the stability or regularity 
of 3d Euler equations. Along these lines, one has the following proposition. 

Proposition 3.1 For 1 < p < oo there is no continuous function r i— > </>(t) 
such that for any smooth solution of the Euler equations the following estimate 

IK", t)\\ w i,p in) < 0(|K-, Q)lw^(Q)) > 

is true. 

Observe that the above statement is not in contradiction with the local stability 
results, which produce local control of higher norm at time t in term of higher 
norm at time as done in (fT7|) according to the formula 



fors>-, \\u(t)\\ HHn) < 



2' 1 v n,n ~ l-Ct|u(0)| H . ( n) ' 

Proof. The proof is done by inspection of a pressureless solution, defined on a 
periodic box (R/Z) 3 of the form 

u(x,t) = (u 1 (x 2 ),0,u 3 {xi - tui(x 2 ),x 2 )), 

which satisfies 

V • u = , d t u + u ■ Vu = . 
Therefore, the initial data satisfies the relation 



IK->°)Ci.p(f2) - J q \d X2 u 1 (x 2 )\ p dx 1 + 

n(\d Xl u 3 (x 1 ,x 2 )\ p dx 1 dx 2 + \d X2 u 3 (x 1 ,x 2 )\ p )dx 1 dx 2 . (34) 

And for t > 

ll' u ("j*)llw , i,j'(n) — / \dx 2 ui(x 2 )\ p dxidx 2 dx 3 + 



1 r l 



(\d Xl u 3 (xi, x 2 )\ p dxidx 2 + \d X2 u 3 (xi, x 2 )\ p )dxidx 2 

o Jo 

+ tP / / \d X2 u 1 {x 2 )\ p \d Xl u 3 (xi,x 2 )\ p dx 1 dx 2 . (35) 
Jo Jo 

Then a convenient choice of u\ and u 3 makes the left-hand side of (f34|) bounded 
and the term 

t p / / \d X2 u 1 (x 2 )\ p \d xl u 3 (x 1 ,x 2 )\ p dx 1 dx 2 , 
Jo Jo 
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on the right-hand side of (|35[) grows to infinity as t — > oo. The proof is then 
completed by a regularization argument. 

Remark 3.1 With smooth initial data, the above construction gives an example 
of global, in time, smooth solution with vorticity growing (here only linearly) for 
t — ► oo. 

As in the case of the Riccati differential inequality y' < Cy 2 , one can obtain 
sufficient conditions for the existence of a smooth solution during a finite interval 
of time (say < t < T). On the other hand, this gives no indication on the 
possible appearance of blow up after such time. Complicated phenomena that 
appear in the fluid, due to strong nonlinearities, may later interact in such a 
way that they balance each other and bring back the fluid to a smooth regime. 
Such phenomena is called singularity depletion. 

An example which seems to illustrate such cancelation has been constructed 
by Hou and Li [35) . It is concerning axi-symmetric solutions of the 3d Eulcr 
equations form rf(z), which obviously possess infinite energy. 

Specifically, let us start with the following system of integro-differential equa- 
tions with solutions that are defined for (z,t) £ (R/Z) x R+ 

tit + 2ipu z — —2vu , Vt + 2ipv z = u 2 — v 2 + c(t) (36) 

t/j z = v, [ v(z,t)dz = 0. (37) 
Jo 

In (|36p , the z independent function c(t) is chosen to enforce the second relation 
of (I37p . which in turn makes the function tj){z,t) 1— periodic in the z direction. 
As a result one has the following: 

Lemma 3.1 For any initial data (u(z, 0), v(z, 0)) <E C m (R/Z), with m > 1, 
the system i36}) and has a unique global, in time, smooth solution. 

Proof. The proof relies on a global a priori estimate. Taking the derivative 
with respect to z variable gives (using the notation (u z ,v z ) — (u',v')) : 

u' t + 2ipu' z - 2ip z u' = -2v'u - 2vu' 
v' t + 2ij)v' z - 2ip z v' = 2uv! - 2vv' . 

Next, one uses the relation ip z — —v, multiplies the first equation by u, multiplies 
the second equation by v, and adds them to obtain 

±(u 2 z +v 2 z ) t +ip(u 2 z +v 2 z ) z = 0. (38) 

The relation (|38p provides a uniform L°° bound on the z— derivatives of u and 
v. A uniform L°° bound for v follows from the Poincare inequality, and finally 
one uses for u the following Gronwall estimate 

Hz, < \\u(z, o)| iooe *IK^o)) 2 +(«(^o)) 2 l^ . (39) 
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Remark 3.2 The global existence for solution of the system h36\) and {Sty , 
with no restriction on the size of the initial data, is a result of delicate bal- 
ance/ cancelation, which depends on the coefficients of the system. Any modifi- 
cation of these coefficients may lead to a blow up in a finite time of the solutions 
to the modified system. On the other hand, the solutions of the system hSb})~ 
may grow exponentially in time. Numerical simulations performed by J35}/ 
indicate that the exponential growth rate in h39jl may get saturated. 

The special structure of the system (j3"o) -([3"T ]) is related to the 3d axi-symmetric 
Euler equations with swirl as follows. Introduce the orthogonal basis 

er = ('-,D), e e = (-^,-,0), e z = (0,0,1), 
r r r r 

and with the solution of the system (|3T>|) and (|37p , construct solutions of the 3d 
(2 + 1/2) Euler equation according to the following proposition. 

Proposition 3.1 Assume that u(z,t) and ip(z,t) are solutions of the systems 
\3b}) and \3ty , then the function 

U(z, t) = — r ^ ' - e r + ru(z, i)e$ + 2rijj(z, t)e z 
oz 

is a smooth solution of the 3d Euler, but of an infinite energy. Moreover, this 
solution is defined for all time and without any smallness assumption on the 
size of the initial data. 

4 Weak limit of solutions of the Navier— Stokes 
Dynamics 

As we have already remarked in the introduction, for both practical problems, 
as well as for mathematical analysis, it is feasible to consider the Euler dynamics 
as the limit of the Navier-Stokes dynamics, when the viscosity tends to zero. 
Therefore, this section is devoted to the analysis of the weak limit, as v — > , 
of Leray-Hopf type solutions of the Navier-Stokes equations in 2d and 3d. We 
will consider only convergence over finite intervals of time < t < T < oo. 
We also recall that v denotes the dimensionless viscosity, i.e., the inverse of the 
Reynolds number. 

4.1 Reynolds stresses tensor and dissipative solutions 

As above, we denote by f2 an open set in H d . For any initial data, u u (x,0) = 
Uo(x) E L 2 (fl), and any given viscosity, v > 0, the pioneer works of Leray [33] 
and Hopf [53j (see also Ladyzhenskaya [31] for a detailed survey), which were 
later generalized by Scheffer [53], and by Caffarelli, Kohn and Nirenberg [5], 
showed the existence of functions u v and p v with the following properties 

u v G L°°((0, T); L 2 (Q)) n L 2 ((0, T); Hq(Q)), for every T e (0, oo). (40) 
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In addition, they satisfy the Navier-Stokes equations 



d t u v + V • (u v <g) u u ) - v/S.u v + \Ip v = , (41) 
V • u v = , u^Oon^, (42) 

in the sense of distributions. Moreover, such solutions satisfy the "pointwisc" 
energy inequality 

-d t \u v (x,t)\ 2 + v\Vu v {x,t)\ 2 + 

\u (x t)\ 2 

V-((u„®u v )(x,t)-vV l - vK 2 71 ) + V-(p u (x > t)u u (x ) t)) <0 (43) 
or in integrated form 

\d t [ \u u (x,t)\ 2 dx + u [ \Vu u (x,t)\ 2 dx < 0. (44) 
^ Jn Jn 

A pair {u u ,p„} which satisfies (|4T)|) . (|42p and (|4"3")) is called a suitable weak so- 
lution of the Navier-Stokes equations, in the sense of Caffarelli-Kohn-Nirenberg. 
If it satisfies, however, the integrated version of the energy inequality (|44|) in- 
stead of the pointwise energy inequality (|43|) it will then be called a Leray-Hopf 
weak solution of the Navier-Stokes equations. 

In two-dimensions (or in any dimension but with stronger hypothesis on 
the smallness of the size of the initial data with respect to the viscosity) these 
solutions are shown to be smooth, unique, and depend continuously on the 
initial data. Furthermore, in this case, one has equality in the relations (I43| 
and ([44]) instead of inequality. 

Therefore, as a result of the above, and in particular the energy inequality 
(|4"4"|) , one concludes that, modulo the extraction of a subsequence, the sequence 
{u u } converges in the weak—* topology of L°°(R^, L 2 (Q j) to a limit u; and 
the sequence {Vp^} converges to a distribution Vp, as v — > 0, for which the 
following holds 

u e L°°(R^,i 2 (0)),V •« = 0inf2, u-n = on dfl, 

u(x 1 t)\ 2 dx + 2v [ [ \Vu\ 2 dxdt< [ \u a (x)\ 2 dx , 
n Jo Jn Jn 

lim (u v ® u v ) = u®u+ lim ((tt„ — u) ® (u v — u)), (45) 

d t u + V • (u ® u) + lim V • [ u - u v ) O (u - u v ) ) + Vp = . (46) 

Observe that the term 

RT(x,t) = lim(u(x,t) - uJx,t)) ® (u(x,t) - u v (x,t)) (47) 

v— >0 

is a positive, symmetric, measure-valued tensor. In analogy with (see below) 
the statistical theory of turbulence, this tensor may carry the name of Reynolds 
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stresses tensor or turbulence tensor. In particular, certain turbulent regions will 
correspond to the support of this tensor. 

This approach leads to the following questions. 

1. What are the basic properties (if any) of the tensor RT(x,t) ? 

2. When does the tensor RT(x,t) identically equal zero? Or, what is equiv- 
alent, when does the limit pair {u,p} satisfies the Euler equations? 

3. When does the energy dissipation 



tend to zero as v — > 0? 

4. Assuming that {u,p} is a solution of the Euler equations is such a solution 
regular enough to imply the conservation energy? 

Hereafter, we will use the following notation for the L 2 — norm 



Remark 4.1 The tensor RT(x,t) is generated by the high frequency spatial 
oscillations of the solution. This feature will be explained in more details in 
section 16'. 11 Therefore, such behavior should be intrinsic and, in particular, 
independent of orthogonal (rotation) change of coordinates. For instance, in the 
2d case, assuming that the function u is regular, the invariance under rotation 
implies the relation 



where a(x,t) and /3(x,t) are some scalar valued (unknown) functions. Thus, 



Of course, this "soft information" does not indicate whether (3{x,i) is zero or 
not. It also does not indicate whether this coefficient is positive, nor how to 
compute it. But this turns out to be the turbulent eddy diffusion coefficient that 
is present in classical engineering turbulence models like the Smagorinsky or the 
ke models (see, e.g., JJ^, J5^, J61^, and JEE/. 

Remark 4.2 Assume that the limit {u, p} is a solution of the Euler equations 
which is regular enough to ensure the conservation of energy, i.e. \u(t)\ 2 = \uo\ 2 ■ 
Then by virtue of the energy relation |^| ) we have 






the equation [£b] ) becomes 
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and by the weak limit relation 

]imMhu v (t)\ 2 > l\u(t)\ 2 -. 

v — *0 z i 

one has that the strong convergence and the relation 

rt 



liminf^ / \Vu„(s)\ 2 ds = 



hold. The following question was then raised by Onsager [59]: "What is the 
minimal regularity needed to be satisfied by the solutions of the 2d or 3d Euler 
equations that would imply conservation of energy?" . The question was pursued 
by several authors up to the contribution of Eynik '251, and Constantin, E and 
Titi \15jj . Basically in 3d it is shown that if u is bounded in £°°(R|, -ff^(O)) , 
with [3 > 1/3 , the energy 

i J \u(x,t)\ 2 dx 

is constant. On the other hand, arguments borrowed from statistical theory 
of turbulence (cf. section \6.2}) . show that the sequence u v will be, in general, 
bounded in L°°(R+, H*(fy) and one should observe that such a statement does 
not contradict the possibility of decay of energy in the limit as v — ► . 

To study the weak limit of Leray-Hopf solutions of the Navier-Stokes dynamics, 
P.L. Lions and R. Di Perna [46j introduced the notion of Dissipative Solution 
of the Euler equations. To motivate this notion, let w(x,t) be a divergence free 
test function, which satisfies w ■ ft = on the boundary <9fi . Let 

E(w) = d t w + P(w ■ Vw) , (49) 

where P is the Leray-Helmholtz projector (see, e.g., [17] ). Then for any smooth, 
divergence free, solution of the Euler equations u(x,t) in Q , which satisfies the 
boundary condition u ■ n = on d£l , one has: 

d t u + V • (u <g> u) + Vp = , 
d t w + V • (w «) w) + Vg = E{w) , 

+ 2(S(w)(u - w), (u - w j) = 2(E{w), u-w), 
where S(w) denotes, as before, the symmetric tensor 

S{w) = i(Vio + (Vw) T ) . 

By integration in time this gives 

\u{t) - w(t)\ 2 < &tf 2 IS(«OMIo.«i»| u (o) - w(0)| 2 
ft 

+2 / e^ 2ls(w){T) ^ dT (E(w)( S ),(u-w)(s))ds. (50) 
Jo 

The above observation leads to the following definition 
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Definition 4.1 A divergence free vector field 

uew~C(H t ;(L 2 (n)) d ), 

which satisfies the boundary condition u ■ n = on d£l, is called a dissipative 
solution of the Euler equations if for any smooth divergence free vector 

field w , with w ■ n — on dtt , the inequality i50\) holds. 

The following statement is easy to verify, but we mention it here for the sake of 
clarity. 

Theorem 4.1 i) Any classical solution u of the Euler equations ill]) is a dis- 
sipative solution. 

ii) Every dissipative solution satisfies the energy inequality relation 

K*)l 2 <K0)l 2 . (5i) 

Hi) The dissipative solutions are "stable with respect to classical solutions" . 
More precisely, if w is a classical solution and u is a dissipative solution of the 
Euler equations, then one has 

\u{t) - w{t)\ 2 < eti 2 l s M( s )IU ds | u (0) - w{0)\ 2 . 

In particular, if there exists a classical solution for specific initial data, then any 
dissipative solution with the same initial data coincides with it. 

iv) In the absence of physical boundaries, i.e. in the case of periodic boundary 
conditions or in the whole space H d ,d = 2,3, any weak limit, as v — > 0, of 
Leray-Hopf solutions of the Navier-Stokes equations is a dissipative solution of 
Euler equations. 

Proof and remarks. The point i) is a direct consequence of the construc- 
tion. To prove ii) we consider w = as a classical solution. As a result, one 
obtains for any dissipative solution, the relation (|51[) . which justifies the name 
dissipative. Furthermore, it shows that the pathological examples constructed 
by Scheffer [64] and Shnirelman |65j are not dissipative solutions of Euler equa- 
tions. 

For the point iii) we use in (|50|) the fact that w being a classical solution 
implies that E(w) = 0. We also observe that this statement is in the spirit of the 
"weak with respect to strong" stability result of Dafermos [20 for hyperbolic 
systems. 

Next, we prove iv) in the absence of physical boundaries. Let u v be a Leray- 
Hopf solution of the Navier-Stokes system, which satisfy an energy inequality 
in (UU), and let w be a classical solution of the Euler equations. By subtracting 
the following two equations from each other 

d t u v + V • (w„ <g) u„) — vAu„ + \lp v = 
dtw + V • (w ® w) — vAro + = -vAw, 
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and taking the L 2 inner product of the difference with (u u — w) one obtains 



We stress that the above step is formal, and only through rigorous arguments 
one can see the reason for obtaining an inequality in (152p , instead of an equality. 
However, this should not be a surprise because we are dealing with Leray-Hopf 
solutions, u v , of the Navier-Stokes system which satisfy an energy inequality in 
([44]) , instead of an equality. 

Now, to conclude our proof we observe that in the absence of physical bound- 
aries one uses the relation 



— v A(w — u v ){x, t) ■ (w — u v ){x, t)dx = v \ |V(u> — u v )(x. t)\ 2 dx (53) 



and the result follows by letting v tend to zero. 

Remark 4.3 The above theorem states in particular, and in the absence of 
physical boundaries, that as long as a smooth solution of the Euler equations 
does exist, it is the limit, as v — > 0, of any sequence of Leray-Hopf solutions of 
the Navier-Stokes equations with the same initial data. In a series of papers, 
starting with Bardos, Golse and Levermore J^, connections between the notion 
of Leray-Hopf solutions for the Navier-Stokes equations and renormalized solu- 
tions of the Boltzmann equations, as defined by P.L. Lions and Di Perna, were 
established. In particular, it was ultimately shown by Golse and Saint Ray- 
mond [32] that, modulo the extraction of a subsequence, and under a convenient 
space time scaling, any sequence of such renormalized solutions of the Boltz- 
mann equations converge (in some weak sense) to a Leray-Hopf solution of the 
Navier-Stokes system. On the other hand, it was shown by Saint Raymond J68f 
that, under a scaling which reinforces the nonlinear effect (corresponding at the 
macroscopic level to a Reynolds number going to infinity), any sequence (mod- 
ulo extraction of a subsequence) of the renormalized solutions of the Boltzmann 
equations converges to a dissipative solution of the Euler equations. Therefore, 
such a sequence of normalized solutions of the Boltzmann equations converges 
to the classical solution of Euler equations, as long as such solution exists. In 
this situation, one should observe that, with the notion of dissipative solutions 
of Euler equations, classical solutions of the Euler equations play a similar role 
for the "Leray-Hopf limit" and the "Boltzmann limit." 

Remark 4.4 There are at least two situations where the notion of dissipative 
solution of Euler equations is not helpful. 

The first situation is concerned with the 2d Euler equations. Let u e (x,t) be 
the sequence of solutions of the 2d Euler equations corresponding to the sequence 
of smooth initial data u e (x, 0). Suppose that the sequence of initial data u e [x, 0) 
converges weakly, but not strongly, in L 2 (fl) to an initial data u(x, 0), as e — » . 




h 2(S(w)(u u - w), {u v - w)) - 2v(A(u v - w), {u v - w)) 

< 2(E(w), u v — w) — (vAw, {u v — w)) . ( 



(52) 
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Then for any smooth, divergence free vector field w, one has, thanks to &50\) . 
the relation 

\u e (t) - w(t)\ 2 < efo 2||SM(s)IUd,| Ue(0) _ w(Q) |2 

+2/ e^ ns{w)(T) ^ dT {E{w),u t -w){s)ds. (54) 



However, with the weak convergence as e — > 0, one has only 
\u(0) - w(0)\ 2 < liminf \u e (0) - w(0)\ 2 , 

and might not hold at the limit, as e — > 0. To illustrate this situation, we 
consider a sequence of oscillating solutions of the 2d Euler equations of the form 

u £ {x,t) = U{x,t,^^-) + 0(e), 
e 

where the map 9 — > U(x, t, 9) is a nontrivial 1— periodic function. In Cheverry 
U3f a specific example was constructed such that 

u = w — lim u t = I U(x, t, 9)d9 
e -° Jo 

is no longer a solution of the Euler equations. The obvious reason for that ( in 
comparison to the notion of dissipative solution), is the fact that 

u(x,o,^l) 

e 

does not converge strongly in L (ft) ■ 

The second situation, which will be discussed at length below, corresponds 
to the weak limit of solutions of the Navier-Stokes equations in a domain with 
physical boundary, subject to the no-slip Dirichlet boundary condition. 

As we have already indicated in Theorem \4-l\ one of the most important 
features of the above definition of dissipative solution of Euler equations is that 
it coincides with the classical solution of Euler equations, when the latter exists. 
This can be accomplished by replacing w in i5(A) with this classical solution of 
Euler equations. Therefore, any procedure for approximating dissipative solu- 
tions of Euler must lead to, in the limit, to inequality i5U\) . Indeed, in the ab- 
sence of physical boundaries, we have been successful in showing, in an almost 
straight forward manner, in Theorem \4.1\ that the Leray-Hopf weak solutions of 
the Navier-Stokes equations converge to dissipative solutions of Euler equations. 
On the other hand, in the presence of physical boundaries, the proof does not 
carry on in a smooth manner because of the boundary effects. Specifically, in 
the case of domains with physical boundaries, inequality i5S\) leads to 

\ ^ Ul '— W \ + (S(w)(u v - w), (u v - w)) + v J \V(w - u u )\ 2 dx 

< (E(w), u v — w) — c(Amj, (u u — w)) + v j d n u v ■ wda . (55) 

Jon 
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The very last term in if 55) J represents the boundary effect. We will discuss below 
the subtleties in handling this term. 



5 No-slip Dirichlet boundary conditions for the 
Navier— Stokes dynamics 

This section is devoted to the very few available results concerning the limit, as 
v — > 0, of solutions of the Navier-Stokes equations in a domain SI c R d , d — 
2,3 with the homogenous (no-slip) Dirichlet boundary condition u v — on 
dil. This boundary condition is not the easiest to deal with, as far as the zero 
viscosity limit is concerned. For instance, the solutions of the of 2d Navier- 
Stokes equations, subject to the boundary conditions u v ■ n — and V A u v =0, 
are much better understood and much easier to analyze mathematically [5] as 
the viscosity v — ► 0. However, the no-slip boundary condition is the one which 
is more suitable to consider physically for the following reasons. 

i) It can be deduced in the smooth (laminar) regime, from the Boltzmann 
kinetic equations when the interaction with the boundary is described by a 
scattering kernel. 

ii) It generates the pathology that is observed in physical experiments, like 
the Von Karman vortex streets. Moreover, one should keep in mind that almost 
all high Reynolds number turbulence experiments involve a physical boundary 
(very often turbulence is generated by a pressure driven flow through a grid!) 

The problem emerges first from the boundary layer. This is because for the 
Navier-Stokes dynamics, the whole velocity field equals zero on the boundary, 
i.e. u v = on dft, while for the Euler dynamics it is only the normal component 
of velocity field is equal to zero on the boundary, i.e. u ■ n = on dft. Therefore, 
in the limit, as the viscosity v — > 0, the tangential component of the velocity 
field of the Navier-Stokes dynamics, u„, generates, by its "jump", a boundary 
layer. Then, unlike the situation with linear singular perturbation problems, 
the nonlinear advection term of the Navier-Stokes equations may propagate 
this instability inside the domain. 

As we have already pointed out, the very last term in (|55[) . i.e. the boundary 
integral term in the case of no-slip boundary condition, 



is possibly responsible for the loss of regularity in the limit as v — > 0. This is 
stated more precisely in the following. 

Proposition 5.1 Let u(x, f) be a solution of the incompressible Euler equations 




(56) 
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ill fl X (0, T] , with the following regularity assumptions. 

S(u) = -{Vu + (Vu) T ) G L 1 ((0, T);i°°(0)), and 

«eI 2 ((0,T);F(!l)), /or S >l. 

Moreover, suppose that the sequence u v , of Leray-Hopf solutions of the Navier- 
Stokes dynamics (with no-slip boundary condition) with the initial data u v (x, 0) = 
u(x,0), satisfies the relation 

lim z/||P«o(V A u,,)\\ ,1 =0, 

where Pqq denotes the projection on the tangent plane to dfl according to the 
formula. 

P 9 n(V A u v ) = V A u v - ((V A u v ) ■ n)n. 
Then, the sequence u v converges to u in C((0,T); L 2 (VL)). 

The proof is a direct consequence of (|55[) and (f56"|) , with w being replaced by 
u. This Proposition can be improved with the following simple and beautiful 
theorem of Kato which takes into account the vorticity production in the bound- 
ary layer {x € f2 \ d(x, dfl) < v\ , where d(x, y) denotes the Euclidean distance 
between the points x and y. 

Theorem 5.1 Let u(x,t) € W h °°((0,T) x Q) be a solution of the Euler dy- 
namics, and let u v be a sequence of Leray-Hopf solutions of the Navier-Stokes 
dynamics with no-slip boundary condition. 

d t Uy — vAu u + V • (u u <g> u y ) + S7p u — , u u (x,t) — on dfl , (57) 

with initial data u u (x, 0) = u{x, 0) . Then, the following facts are equivalent. 

(i) lim v [ ( (V Au v ) ■ (ft Au)dadt = (58) 

Jo Jon 

(ii) U v (t) — > u(t) in L 2 (Q) uniformly int G [0,T] (59) 

(iii) u v (i) — > u(i) weakly in L 2 (il) for each t £ [0,T] (60) 

(w) Yixav { [ \Vu u (x,t)\ 2 dxdt = (61) 
"-* Jo Jn 

(v) Yanv j / | Vu„(x, t)\ 2 dxdt = . (62) 



i/->0 



On{d(a;,90)<iy} 



Sketch of the proof. The statement (|59")l is deduced from (1551) by replacing 
w by u in (|55|) and ([56")) . No proof is needed to deduce (f6"0"| from ([59")) , or (|6"Tj) 
from (|6"0)) . 
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Next, one recalls the energy inequality (j4"4")l , satisfied by the Leray-Hopf 
solutions of the the Navier-Stokes dynamics 



If, , _.,o . /•' /• _ .„ . . 1 



2 



\u u {x,T)Ydx + v / / |Vu„(a;,t)rda;dt < - / |«(ar, 0)pda; . (63) 
'o Jn 



By virtue of the weak convergence, as stated in (|60| . and the fact u is a smooth 
solution of the Euler dynamics, one has 

lim i / |«„(a:, T)| 2 dac > - / | lim uJx, T)\ 2 dx 
u^o 2 J n 2 J n u^o 

= \( \u(x,T)\ 2 dx=\ [ \u(x,0)\ 2 dx. (64) 
1 Jn 1 Jn 

Together with ([53]) this shows that implies |[6"T]). 

The most subtle part in the proof of this theorem is the fact that ([62]) implies 
(|58[) . The first step is the construction of a divergence free function v v (x, t) with 
support in the region {x G f2|d(:r, 9fi) < v\ x [0, T), which coincides with u on 
d£l x [0,T], and which satisfies (with K being a constant that is independent of 
v) the following estimates 

IKIU°°(nx(Q,T)) + \\d(x, 9fl)Vu v |io<.(nx(Q,T)) < # 5 (65) 
||(d(x,9r!)) 2 V^|| ioo(0x(0 , T)) < JiV, (66) 

IKli°°((0,T);i 2 (n)) + l|5t^|L~((0 : T);L2(O)) < , (67) 

||Vu„||ioo(( 0i r);i2(n)) < , (68) 

||Vu v li«( n>< (o,r)) < • (69) 

Then we multiply the Navier-Stokes equations by v K and integrate to obtain 

pT p pT p 

—v \ I (V A u„) • (n A u)dadt = —v \ I ■ v^dadt 

Jo Jan Jo Jan on 

= —v I j Au v ■ v v dxdt — v \ j (Vu v : \7v u )dxdt (70) 
Jo Jn Jo Jn 

T pT pT 

{d t u v ,v v )dt- I (V • {u v ®u v ),v v )dt-v I {Vu„,Vv v )dt . (71) 



Eventually, by using the above estimates and (|62p , one can show that 

l^ Q (j {(d t u v ,v u ) + (V-(u u ®u v ),v v ) + v(Vu l/ ,Vv v ))dtj =0, (72) 
which completes the proof. 

Remark 5.1 In Constantin and Wu J18}/ the authors study the rate of con- 
vergence of solutions of the 2d Navier-Stokes equations to the solutions of the 
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Euler equations in the absence of physical boundaries, for finite intervals of time. 
Their main observation is that while the rate of convergence, in the L 2 -norm, 
for smooth initial data is of the order Oiy) it is of the order of 0{^/v) for less 
smooth initial. The order of convergence 0{^fv) is, for instance, attained for 
the vortex patch data with smooth boundary. In this case the fluid develops an 
internal boundary layer which is responsible for this reduction in the order of 
convergence. 

In the 2d case and for initial data of finite W 1,p norm, with p > 1, (and 
also for initial data with vorticity being a finite measure with a "simple" change 
of sign J21f . \49p one can prove the existence of "weak solutions of the Eu- 
ler dynamics. " In the absence of physical boundaries, such solutions are limit 
points of a family of (uniquely determined) Leray-Hopf solutions of the Navier- 
Stokes dynamics. However, these weak solutions of the, Euler equations are not 
uniquely determined, and the issue of the conservation of energy for these weak 
solutions is, to the best of our knowledge, completely open. 

On the other hand, in a domain with physical boundary and with smooth 
initial data, Theorem \5.1\ shows a clear cut difference between the following two 
situations (the same remark being valid locally in time for the 3d problems). 

i) The mean rate of dissipation of energy 

e= - jf J \Vu v (x,t)\ 2 dxdt 

goes to zero as v — > 0, and the sequence u v of Leray-Hopf solutions converges 
strongly to the regular solution u of the Euler dynamics. 

ii) The mean rate of dissipation of energy does not go to zero as v — ► 
(modulo the extraction of a subsequence) , so the corresponding weak limit u of 
u v does not conserve energy, i.e. 

\\u(t)\ 2 < i| U (0)| 2 , for some t€(0,T), 

and one of the following two scenarios may occur at the limit. 

a) At the limit one obtains a weak solution (not a strong solution) of the 
Euler dynamics that exhibits energy decay. Such a scenario is compatible with 
a uniform estimate for the Fourier spectra 

E v {k,t) = \u v {k,t)\ 2 \k\ d - 1 , fi„(fc > t) = _J_ / e - ikx u v (x,t)dx , 

which may satisfy a uniform, in v, estimate of the following type 

E v (k,t)<C\k\-P , (73) 

provided f3 < 5/3. Otherwise, this would be in contradiction with the results of 
Onsager i59\j . Eyink \25^ . and Constantin, E and Titi \151 - 

b) No estimate of the type J7ff[ ) is uniformly (in the viscosity) true, and 
the limit is not even a solution of the Euler dynamics, rather a solution of a 
modified system of equations with a term related to turbulence modeling - an 
"eddy-viscosity" like term. 
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6 Deterministic and statistical spectra of tur- 
bulence 



6.1 Deterministic spectra and Wigner transform 

The purpose of this section is the introduction of Wigner measures for the 
analysis (in dimensions d = 2,3) of the Reynolds stresses tensor 

RT(u v )(x, t) = lim ((it„ - u) <g> (u v - u)) , (74) 

v— >0 

which appears in the weak limit process of solutions of the Navier-Stokes equa- 
tions it„, as v — > 0, cf. ([431) . (Notice that RT(u v )(x,t) is independent of the 
viscosity v, but it depends on the sequence {u v }.) This point of view will be 
compared below (cf. section 16. 2[) to ideas emerging from statistical theory of 
turbulence. 

Let {u„,pv} be a sequence of Leray-Hopf solutions of the Navier-Stokes 
equations, subject to no-slip Dirichlet boundary condition in a domain f2 (with 
a physical boundary). Thanks to the energy inequality 1|44|) (possibly equality 
in some cases) 



±\u v {- : t)\ 2 + v \Vu v (-,t)\ 2 dt< ^K-,0)! 2 , (75) 

the sequence {«„} converges (modulo the extraction of a subsequence), as v — > 0, 
in the weak—* topology of the Banach space L°°((0, T); L 2 (il)), to a divergence 
free vector field u, and the sequence of distributions {Vp v } converges to a 
distribution Vp. Moreover, the pair {u, p} satisfies the system of equations 

V • u = in , u ■ n = on dft , (76) 
<9 t u + V • (u ® u) + V • iJT^) + Vp = in O , (77) 

(cf. (|46| V Being concerned with the behavior of the solution inside the domain, 
we consider an arbitrary open subset Q' , whose closure is a compact subset of 
fl, i.e. W CC ft. Assuming the weak—* limit function u belongs to the space 
^((O^T);!! 1 ^)) we introduce the function 

v v = a{x){u u — u) , 

with a(x) e X>(f2) , a(x) = 1 for all x £ 57' . As a result of ((75)) . and the above 
assumptions, the sequence v v satisfies the uniform estimate 

v [ f \Vv„\ 2 dx <C. (78) 
Jo Jn 

Consequently, the sequence v v is (in the sense of Gerard, Mauser, Markowich 
and Poupaud |30j ) \/v— oscillating. Accordingly, we introduce the deterministic 
correlation spectra, or Wigner transform, at the scale : 

RT^~ u )(x,t,k) = --^ 3 [ e ik -v(v u (x-y^-y)®v v (x+y^y))dy. 

(Stt)" 2 2 
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By means of the inverse Fourier transform, one has 



v v (x,t) <g> v v {x,t) = / RT(v v )(x,t,k)dk. (79) 



The tensor RT(v u )(x,t,k) is the main object of section 1 of [30]. Modulo the 
extraction of a subsequence, the tensor RT{y v )(x, t, k) converges weakly, as 
v — > 0, to a nonnegative symmetric matrix- valued measure RT(x,t,dk), which 
is called a Wigner measure or Wigner spectra. Moreover, inside the open subset 
Q' the weak limit u is a solution of the equation 

d t u + V • (u <x> u) + V • / (x, t, djfe) + Vp = . 

The Wigner spectra has the following properties 

i) It is defined by a two points correlation formula 

ii) It is an (x,t) locally dependent object. Specifically, for any <fi £ T>(Q,) one 
has: 

lim (j^y J Rd e ik - v ((^)(x - ^-V) ® (<Pv„)(x + ^-y))dy^ 

= \(p(-)\ 2 RT(x,t,dk) (80) 

Therefore, the construction of RT{x, t, dk) is independent of the choice of the 
pair a(x) and the open subset 0'. 

iii) It is a criteria for turbulence: Points (x, t) around which the sequence w„ 
remains smooth and converges locally strongly to u, as v — > 0, are characterized 
by the relation 

Trace(fif (x, t, dk)) = . (81) 

iii) It is a microlocal object. In fact, it depends only on the behavior of the 
Fourier spectra of the sequence <j>{x)v v (or in fact (j){x)u v ) in the frequency 
band 

A< \k\ < 4=- 



Proposition 6.1 For any pair of strictly positive constants (A, B) and any test 
functions {tj}{k),^{x),6{t)) E C^(Rf) x Cg°(R£) x C^(R+) one has 

V>(fc)|^(x)| 2 0(t)Trace(i?T(^))(a;, t, dk)dxdt 

lim / 6(t) f {^{^k)X^,{^))dkdt . (82) 



A<|fc|<^= 
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The only difference between this presentation and what can be found in [3D] 
comes from the fact that the weak limit u has been subtracted from the sequence 
u v . Otherwise the formula (|T5|) together with the energy estimate are the first 
statement of the Proposition 1.7 of [3D]; while the formula (|82p is deduced from 
(1.32) in [3DJ by observing that the weak convergence of u v — u to implies that 

lim ( / 0(t) f &(Vvk)J^,J4w^)dkdt) =0. 



i/->0 



k\<A J 



6.2 Energy spectrum in statistical theory of turbulence 

The Wigner spectra studied in the previous section turns out to be the deter- 
ministic version of the "turbulent spectra" , which is a classical concept in the 
statistical theory of turbulence. The two points of view can be connected with 
the introduction of homogenous random variables. Let (9JI, 5, dm) be an under- 
lying probability space. A random variable u(x, fj,) is said to be homogenous if 
for any function F the expectation of F(u(x, //)), namely, 

(F{u(x,-))) = [ F{u{x,n))dm{n) 

is independent of x, that is 

V x ((F(u(x, ■)))=(). 

In particular, if u(x,fi) is a homogeneous random vector-valued function, one 
has 

(u(x + r, •) ® u(x, ■)) = (u(x + -, •) <g> u(x - -, •)) , (83) 
which leads to the following proposition. 

Proposition 6.1 Let u(x,/j,) be a homogenous random variable and denote by 
u(k,[i) its Fourier transform. Then one has: 

r , , r 



(fc, .) <8 .)) = __ / e - lk ' r {u[x + -,-)®u(x- -,-))dr. (84) 



(2n) d J Rd " x T ' v 2' 



Proof. The proof will be given for a homogenous random variable which is 
periodic with respect to the variable x, with basic periodic box of size 2nL. The 
formula (|84|) is then deduced by letting L go to infinity. From the Fourier series 
decomposition in (R/27rLZ) d , one has 



u(k, fi) (g) u(k, /i) = 

TTTxId / {u{y^)e' 1 ^ ®u{y'^)e ih ^)dy'dy^ 

\l-KL) J(R/2 7 rI,Z)< i J(TL/2nLZ) d 

! ' ~~ 1 (u(y, O u(y + r, n))dydr . 



(2-KL) 2d J( R / 2 TrLZ) d i(R/27rLZ 
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Averaging with respect to the probability measure dm, using the homogeneity 
of the random variable u(x, [i) and integrating with respect to dy, gives us the 
following 



(fi(fc,.)®fi(fc,.)> = — —r-f / e 1 l ( Uv ( y+ -,.)<8>u„(y- -,.))dr. (85) 

This concludes our proof. 

Next, assuming that in addition to homogeneity, the expectation of the two 
points correlation tensor, (u(x + r, •) ® u(x,-)), is isotropic (i.e., it does not 
depend on the direction of the vector r, but only on its length) one obtains the 
following formula 

(u(k,-)®u(k,-)) = 1 f e- lh ^{u u {y+ r - r )®u v {y- r -,-))dr 
(2nL) a J( R/2TrLZ)t i 2 2 

E(\k\) k®k 
S d -!\k\ d -^ \k\ 2 h 

with Si = 2n , 5*2 = 4-7T, which defines the turbulent spectra 2?(|fc|) ■ 

The notion of homogeneity implies that solutions of the Navier-Stokes equa- 
tions satisfy a local version of the energy balance, often called the Karman- 
Howarth relation cf. (|86|) . Specifically, let {u Ul p u } be solutions of the forced 
Navier-Stokes equations in f2 (subject to either no-slip Dirichlet boundary con- 
dition, in the presence of physical boundary, or in the whole space or in a 
periodic box) 

d t u v + V • (u v (g) u v ) - vAu v + Vp u = f . 

Here u v and p v are random variables which depend on (x,t). We will drop, 
below, the explicit dependence on \i when this does not cause any confusion. 

Multiplying the Navier-Stokes equations by u v {x,t,pi) and assuming that 
one has the following equality 

^-d t \u u (x,t, y)\ 2 - \7 X ■ {{v\7 x u v - p v I)u v )(x,t, fi) + 
u\V x u u (x, t, /i)| 2 = f(x, t) ■ u v (x, t, /j.) , 

we observe, that in the 2d case, the above relation is a proven fact. However, 
in the 3d case, the class of suitable solutions of the Navier-Stokes equations, in 
the sense of Caffarelli-Kohn-Nirenberg, are known to satisfy a weaker form of 
the above relation, involving an inequality instead of equality (cf. (|43[) ). 
Thanks to the homogeneity assumption, the quantity 

(((fV x u„ -p u I)u v )(x,t,-)) 

does not depend on x, and therefore 

(V x • {{{vV x u v - p v I)u v ))(x,t, ■)) = V x ■ {{{vVxUv - p v I)u v )(x,t, ■)) = 0. 
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Thus, the averaged pointwise energy relation 



~dt(\u v (x, t, -)| 2 > + u{\Vu v {x, t, -)| 2 ) = (f(x, t) ■ u v (x, t, •)) (86) 

is obtained. The above is often called the Karman-Howarth relation, and it 
implies that the quantities 



(\u v (x,t,-)\ 2 ) and e{v) = j / k{\T>i„[.v. s..)\-},ls 



t 

2\ 





are uniformly bounded in time, under reasonable assumptions on the forcing 
term. 

Finally, assume also that the random process is stationary in time. Then, 
the equation (]86|) gives an a priori estimate for the mean rate of dissipation of 
energy 

e(v)=v{\\7u v (x,t,-)\ 2 ), 

which is independent of (x,t). Now, with the forcing term / acting only on 
low Fourier modes, one may assume the existence of a region, called the inertial 
range, where the turbulence spectra -E(|fc|) behaves according to a universal 
power law. The name inertial refers to the fact that in this range of wave 
numbers, the energy cascades from low modes to high modes with no leakage of 
energy. That is, there is no viscous effects in this range and only the inertial 
term (u- V)it is active. Combining the above hypotheses: homogeneity, isotropy 
and stationarity of the random process, together with the existence of an inertial 
range of size 

A < \k\ < Bv'i , 

where the spectra behaves according to a power law, one obtains finally, by a 
dimensional analysis and in the three-dimensional case, the "Kolmogorov law" : 

E(\k\)~e(v)l\k\-i . (87) 

It is important to keep in mind the fact that this derivation is based on the 
analysis of a random family of solutions. Therefore, the formula (|87|) combined 
with the formula (|85|) implies that in the average the solutions have a spectra 
which behaves in the turbulent regime according to the prescription. 

1 f -k-r r r 

{u v (k,t,-)®u v {k,t,-)) = — \ e- t —(u u {y+-,t,-)<S>u„(y--,t, 

^ e(f)J|fe| _ t k <g> k 
~ Airl^- 1 1 ~ ~\kY 

Remark 6.1 The main difficulty in the full justification of the above derivation 
is the construction of a probability measure, dm, on the ensemble of solutions of 
the Navier-Stokes equations that would satisfy the hypotheses of homogeneity, 
isotropy and stationarity. In particular, the construction of such measure should 



28 



be uniformly valid when the viscosity v tends to . See, for instance, the books 
of Vishik and Fursikov and Foias et al. f28jj for further study and references 
regarding this challenging problem. 

The next difficulty (which is a controversial subject) is the justification for 
the spectra of an inertial range with a power law. In Foias et al. J28]j (see also 
Foias 1271) it was established, for example, the existence of an inertial range 
of wave numbers where one has a forward energy cascade. However, we are 
unaware of a rigorous justification for a power law in this inertial range. 

Nevertheless, the construction of a power law of the spectra is often used as a 
benchmark for validation of numerical computations and experiments. Since, in 
general, one would have only one run of an experiment, or one run of simulation, 
a Birkhoff theorem, which corresponds to assuming an ergodicity hypothesis, is 
then used. This allows for the replacement of the ensemble average by a time 
average. For instance, one may assume, in presence of forcing term, in addition 
to stationarity that for almost all solutions, i.e. almost every /i, one has 

1 f T r r 

lim — / u u (y+-,t,(jL)®u 1 ,(y- -,t,p)dt = 

1 —*oo 1 Jq Z Z 

r r 

which would give the following relation, for almost every solution u v , 
1 ' T 



lim — / u v (k,t, fi) ® u u (k,t, fj,)dt 

T — >oo _Z 







7^ I e ik - r {u v (y + T -,t, •) <g> u v {y - T -,t, -))dr 
\zit) a jRd 2 2 

e{v)i \k\~3 k®k^ 



= (u(k, t, •) (8) u(k, t, ■)) ~ \^ d U {I ~ -j£p-)> • 

6.3 Comparison between deterministic and statistical spec- 
tra. 

The deterministic point of view considers families of solutions u v of the Navicr 
Stokes dynamics, with viscosity v > 0, and interprets the notion of turbulence 
in terms of the weak limit behavior (the asymptotic behavior of such sequence 
as v — > 0) with the Wigner spectra: 

RT(x,t,dk) = 

lim n?Ad / e ik -y((u„-u){x- ^y,t)®{u v -u){x + ^-y,t))dy. 

As we have already observed earlier, this is a local object (it takes into account 
the (x,t) dependence). Moreover, one could define the support of turbulence, 
for such family of solutions, as the support of the measure Tracei?T(x, t, dk). 
Of course, determining such a support is extremely hard and is a configuration 
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dependent problem. This is perfectly described in the sentences of Leonardo da 
Vinci, who is very often quoted, and in particular in page 112 of [29j : 



doue la turbolenza dellacqua sigenera 

done la turbolenza dellacqua simantiene plugho 

doue la turbolenza dellacqua siposa. 

Up to this point nothing much can be said without extra hypotheses, except 
that the formula (|82|) indicates the existence of an essential, if not an "inertial" , 
range 

A<|fc|<^. (88) 

On the other hand, the statistical theory of turbulence starts from the hy- 
potheses, that seem difficult to formulate in a rigorous mathematical setting, 
concerning the existence of statistics (a probability measure) with respect to 
which the two points correlations for any family of solutions, u v , of the Navier- 
Stokes equations are homogeneous and isotropic. Under these hypotheses, one 
proves properties on the decay of the spectra of turbulence. Moreover, with all 
these assumptions one obtains, by simple dimension analysis, for averages of 
solutions with respect to the probability measure dm(/J,) , the following formula 
in the inertial range: 

(MM, -)l 2 > = (^)i^r * ■ (89) 

Finally, by assuming and using the stationarity (in time) of these two points 
correlations, and by applying, sometimes, the Birkhoff ergodic theorem, one 
should be able to obtain, for almost every solution, the formula 

1 " T 



lim — / (u„(k,t, ® u„(k,t, fi))dt 







(e(i/))3|fc|~3 k®k, 



{u v {k, t, •) ® a„(fc, t, •)) ^ 47r | fc |rf-i ( J " -ftp-)) 



Remark 6.2 Some further connections between these two aspects of spectra 
may be considered. 

i) Assuming that near a point (xq, to) the Wigner spectra is isotropic. Define 
the local mean dissipation rate of energy as 

t)\ 2 dxdt, (90) 

Jo Jn 

with (f> being a localized function about {xa,to). Then one can prove that, for \k\ 
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in the range given by \88\l . 



(27r) d 



1 





Hx a ,t )i v )\ k \ 'hi 



k®k 




as v tends to zero. 

ii) Give sufficient conditions that will make the Wigner spectra isotropic. 
This is in agreement with the fact that this spectra involves only a small scale 
phenomena; thus, it is a reasonable hypothesis. However, the example con- 
structed by Cheverry U3]j shows that this has no chance of always being true. 

Hi) Another approach for establishing the existence of an inertial range for 
forward energy cascade in id, and forward enstrophy cascade in 2d, is presented 
in J28f (see also references therein). This approach is based on the statistical 
stationary solutions of the Navier-Stokes equations. These are time indepen- 
dent probability measures which are invariant under the solution operator of the 
Navier-Stokes equations. Furthermore, in Foias \27^ some semi-rigorous argu- 
ments are presented to justify the Kolmogorov power law of the energy spectrum. 

7 Prandtl and Kelvin— Helmholtz problem 

In this section it is assumed that the sequence of solutions {u v } of the Navier- 
Stokes equations with no-slip Dirichlct boundary condition (in the presence of 
physical boundary) converges to the solution of the Euler equations. According 
to Theorem 15. II of Kato, in this situation one has 



However, since the tangential velocity of the solution of the Euler equations 
is not zero on the boundary, as v — ► 0, a boundary layer is going to appear. On 
the one hand, the scaling of the boundary layer has to be compatible with the 
hypothesis (j9"Tj) ; and on the other hand the equations that model the behavior 
in this boundary layer have to reflect the fact that the problem is very unstable. 
This is because the instabilities (and possible singularities) that occur near 
the boundary may not remain confined near the boundary, and will in fact 
propagate inside the domain by the nonlinear advection term of the Navier- 
Stokes equations. These considerations explain why the Prandtl equations (PE) 
of the boundary layer are complicated. 

There are good reasons to compare the Prandtl equations with the Kelvin- 
Helmholtz problem (KH): 

1. Even though some essential issues remain unsolved for KH, it is much 
better understood from the mathematical point of view than the PE problem. 
However, the two problems share similar properties such as instabilities and 
appearance of singularities. 




(91) 
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2. At the level of modeling, in particular for the problem concerning the 
wake behind an air plane and the vortices generated by the tip of the wings, it 
is not clear if turbulence should be described by singularities in KH or PE (or 
both)! 

For the sake of simplicity, these problems are considered in the 2d case, and 
for the PE in the half space > . 

7.1 The Prandtl Boundary Layer 

One starts with the 2d Navier-Stokes equations in the half plane X2 > 0, with 
the no-slip boundary condition u v {xi, 0, t) = 



d t v!i - i/Av? + u\d Xl u\ + u v 2 d X2 u\ + d xlP v = , (92) 

dtv,2 — vIS.u v 2 + Uid Xl U2 + u v 2 d X2 u2 + d X2 p u = , (93) 

d Xl u v i + d*X 2 = (94) 

<(xi,0,i) = i^(xi,0,t) = on xi e R, (95) 



and assumes that inside the domain (away from the boundary) the vector field 
u I/ (xi, X2-, t) converges to the solution Ui n t(£i, X2, t) of the Euler equations with 
the same initial data. The tangential component of this solution on the bound- 
ary, X2 = 0, and of the pressure are denoted by 

U{ Xl ,t) =uf\xi, 0,t), P(x u t) =p(x u 0,t). 

Then one introduces the scale e = \/v. Taking into account that the normal 
component of the velocity remains on the boundary, one uses the following 
ansatz, which corresponds to a boundary layer in a parabolic PDE problem 

Inserting the right hand side of (|96| into the Navier-Stokes equations, and 
returning to the notation {x\,x?j) for the variables 

X2 

Ai = xi,X 2 = — , 
e 

and letting e go to zero, one obtains formally the equations: 

Ui(xi,0,t) + Ui(xi,0,t) =0, (97) 

d X2 p(xi,x 2 ) = 0^ p{xi,x 2 ,t) = P{xi,t) (98) 

d t u\ - dl 2 ui + u x d Xl ui + u 2 d X2 ui = d Xl P(xi,t) , (99) 

d Xl ui + d X2 U2 — , ui(xi, 0) = t2 2 (a;i, 0) = for x% S R , (100) 

lim Ui(x%, X2) = lim U2(%i, ^2) = . (101) 
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Remark 7.1 As an indication of the validity of the Prandtl equations we oh 
serve that A96\) is consistent with Kato Theorem \5.1[ Specifically, thanks to {Pi 
one has 

T , 

\Ju u {x,s)\ 2 dxds < Cy/v. 

an{d(x,dsi)<cv} 

Remark 7.2 The following example, constructed by Grenier J33f , shows that 
the Prandtl expansion cannot always be valid. In the case when the solutions 
are considered in the domain 

(R 9!1 /Z)xH+ . 

Grenier starts with a solutions u" e j of the pressureless Navier-Stokes equations 
given by 

Kef = Ke/(WVH°) 

d t U re f - dyyUref = , 

where Y — y / \Jv . Using a convenient and explicit choice of the function u re f , 
with some sharp results on instabilities, a solution of the Euler equations of the 
form 

u = u ref + Sv + 0{S 2 e 2Xt ) for < t < ' 



A log 5 

is constructed. It is then shown that the vorticity generated by the boundary for 
the solution Navier-Stokes equations (with the same initial data) is too strong 
to allow for the convergence of the Prandtl expansion. One should observe, 
however, that once again this is an example which involves solutions with infinite 
energy. It would be interesting to see if such an example could be modified to 
belong to the class of finite energy solutions; and then to analyze how the modified 
finite energy solution might violate the Kato criteria mentioned in Theorem \5.1\ 

It is important to observe that in their mathematical properties the PE exhibit 
the pathology of the situation that they are trying to model. First one can prove 
the following proposition. 

Proposition 7.1 Let T > be a finite positive time, and let (U(x,t), P(x, t)) £ 
C + "(Rt x (R Xl x R+ )) be a smooth solution of the 2d Euler equations satis- 
fying, at time t = 0, the compatibility condition Ui(xi,0,t) = [^(asi, 0, i) = 
(notice that only the boundary condition U2(xi,0,t) = is preserved by the 
Euler dynamics). Then the following statements are equivalent: 

i) With initial data u(x, 0) = , the boundary condition Ui(xi, 0, t) = Ui{x\,Q, t) 
in \9T) , the right-hand side in Wtfy given by P(xi,t) = P(xi,0, t) , and the 
Prandtl equations have a smooth solution u(x,t), for < t < T . 

ii) The solution u u (x,t), of the Navier-Stokes equations with initial data 
u v {x,Q) = U(x,0) and with no-slip boundary condition at the boundary xi = 0, 
converges in C 2+a to the solution of the Euler equations, as v — > 0. 
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The fact that statements i) and ii) may be violated for some t is related to the 
appearance of a detachment zone, and the generation of turbulence. This is 
well illustrated in the analysis of the Prandtl equations written in the following 
simplified form 

d t ui - dl 2 ui + uid Xl ui + u 2 d X2 ui = , d Xl ui + d X2 u 2 = (102) 
ui(x 1( 0) = u a (xi,0) = for x l G R (103) 
lim ui(xi,x 2 ) = 0, (104) 

X 2 — >00 

Ui{xi,X2,Q) = uo(xi,x 2 ) ■ (105) 

Regularity in the absence of detachment corresponds to a theorem of Olcinik 
[55]. She proved that global smooth solutions of the above system do exist 
provided the initial profile is monotonic, i.e. for any initial profile satisfying 

u(x,0) = (ux(xi,X2),0) , d X2 ui{xi,x 2 ) ^ 0. 

On the other hand, initial conditions with "recirculation properties" leading to 
a finite time blow up have been constructed by E and Engquist [33J and [2"1] . 
An interesting aspect of these examples is that the blow up generally does not 
occur on the boundary, but rather inside the domain. 

The above pathology appears in the fact that the PE is highly unstable. 
This comes from the determination u 2 in term of u\ by the equation 

d Xl iti + d X2 u 2 = 0. 

Therefore, it is only with analytic initial data (in fact analytic with respect to 
the tangential variable is enough) that one can obtain (using an abstract version 
of the Cauchy-Kowalewskya theorem) the existence of a smooth solution of the 
Prandtl equation for a finite time and the convergence to the solution of the 
Euler equations during this same time (Asano pQ, Caflisch-Sammartino [10] . 
and Cannone-Lombardo-Sammartino |12j.) 



7.2 The Kelvin-Helmholtz problem 

The Kelvin-Helmholtz (KH) problem concerns the evolution of a solution of the 
2d Euler equations 

d t u + V • («® u) + Vp = 0, V-u = (106) 

with initial vorticity w(x,0) being a measure concentrated on a curve T(0) . 

This is already simpler than the PE because the pathology, if any, should in 
principle be concentrated on a curve. Furthermore, the dynamics in this case 
inherits the general properties of the 2d dynamics. In particular, it will obey 
the equation 

<9 t (V A«) + (ii' V)(V A u) = , 

for the conservation (for smooth solutions) of the density of the vorticity. There- 
fore, one can guarantee the existence of a weak solution when the initial vorticity, 
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w(x,0) , is a Radon measure. This was done first by Delort, assuming that the 
initial measure has a distinguished sign [21] . Then the result was generalized to 
situations where the change of sign was simple enough [35]. However, this re- 
markable positive result is impaired by the non-uniqueness result of Shnirelman 

For smooth solutions of the KH, i.e. the ones with vorticity ui - a bounded 
Radon measure with support contained in curve St = {r(A,t), A € R} - the 
velocity field is given, for x ^ E t , by the so called Biot-Savart law 

1 f x — r' 

U'{x,t) = —R^ / | m w(r',t)ds' 

2ir 2 J \x — r'\ z 

1 „ f x-r(X',t) , x ds(X',t) , w . 

: = J w^kw^^^^- (107) 

Here r(A,t), for A G R, is a parametrization of the curve S t , and s(A, t) = 
|r(A, t)| is the corresponding arc length. denotes the -|— counterclockwise 
rotation matrix 

' -1 
1 



Ri 



Furthermore, as x approaches the curve S t the velocity field u admits the 
two-sided limits u± . By virtue of the incompressibility condition one has the 
continuity condition for the normal component of the velocity field, i.e. 



w_ ■ n = Uj 



where hereafter r and n will denote the unit tangent and unit normal vectors 
to the curve E t , respectively. In addition, the average 



. . u + + U- 
(u) = 



2 

is given by the principal value of the singular integral appearing in ()107p 



' = " = p- v / ^p^'' *W ■ (los) 



Using the calculus of distributions one can show that, as long as the curve £ t is 
smooth, the velocity field u, defined above, being a weak solution of the Euler 
equations 

d t u + V • (u <g> u) + X7p = , V • u = , 

is equivalent to the vorticity density u> and the curve S/ satisfying the coupled 
system of equations 

ut - d s (u(d t r - v) ■ = , (109) 
(r t -v)-n = Q, (110) 
"(^) = ^!P- V / T^2"(r',t)ds'. (Ill) 
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The equations (|109[) . (|1 10[) and pil[) do not completely determine r(X,t) . This 
is due to the freedom in the choice of the parametrization of the curve E t . 
Assuming that w^O one introduces a new parametrization X(t, s) which reduces 
the problem to the equation 

or with the introduction of the complex variable z = r± +ir^, where r = (rj., r%), 
one obtains the Birkhoff-Rott equation 

w-aW ' ^«)%.o - (113) 

Remark 7.3 TTie following are certain mathematical similarities of the KH 
problem with the PE: 

1. As for the PE one has for the evolution equation \113\) a local, in time, 
existence and uniqueness result in the class of analytic initial data. This is 
done by implementing a version of the Cauchy-Kowalewskya theorem (Bardos, 
Frisch, Sulem and Sulem \61^ .) 

2. As for the PE one can construct solutions that blow up in finite time. 

3. One observes that the singular behavior in the experiments and numerical 
simulations with the KH problem is very similar to the one that is generated by 
the no-slip boundary condition when the viscosity is approaching zero. 

The best way to understand the structure of the KH is to use the fact that 
the Euler equations are invariant under both space and time translations, and 
under space rotations, and to consider a weak solution of the 2c? Euler dynamics 
either in the whole plane R 2 , satisfying 

ueC((-r,T);L 2 (R 2 )), T>0 

or subject to periodic boundary conditions satisfying 

ueC((-T,T);L 2 ((R/Z) 2 )). 

Assuming that in a small neighborhood U of the point (t = 0, z = 0) the vorticity 
is concentrated on a smooth curve in the complex plane which takes the form: 

z(\,t) = (at + {3(\+ef(\,t)), /(0,0) = V/(0,0) = 0. (114) 

Then using the relations V-ii = 0,VAu = w and the Biot-Savart law, one 
obtains 

e\(3\ 2 dj(\,t) = 

^"P.V.f dA f ,-„ J7TTT, +E(z(X,t)) (115) 

2™ J M t,y)eu} (A-A0(l-6 /(A 'ty ,0 ) 
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where here, and in the sequel, E(z) denotes the "remainder" , which is analytic 
with respect to z. Next, use the expansion 



2tt 



1 p - v / { x-x>)(i-%tr^ dX ' = (116) 

f(X,t)-f(X',t)^, / , (/(A,t)-/(A / ,t))" dA/ 



[ f(\,t)-f(\',t) ^, ^ e" /■ 



(A - A')(" +1 ) 

and implement the following formulas concerning the Hilbert transform 

to deduce, from (|115[) and (|116p . that the real and imaginary parts of f(X,t) — 
X(X, t) + iY(X, t) are local solutions of the system 



1 

21 



d * X = ^2 \ D * \ Y + ei? i ( X > Y ) + ( X > Y ) (! 19 ) 



or 



d t Y=-^\D x \X + eR 2 (X,Y) + E 2 (X,Y) (120) 



(dt + ^ d D Y = <d t Ri{X,Y) - ^—\D X \R 2 {X,Y)) 
+d t E 1 (X,Y)-^\D x \E 2 (X,Y) (121) 
( 9 t + 7^u d D Y = <d t R 2 {X,Y) - J—^R^Y)) 



m 4 

i 
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+d t E 2 {X,Y)- 1 ^—\D x \E l {X,Y). (122) 



In ([T2T1 ) and (TT22")) the terms 

dtE^Y) - * |£> A |ifc(*,y) and d t E 2 (X,Y) - ■^^\D\\Ei(X, Y) 

are the first order derivatives of analytic functions with respect to (X, Y) while 
the terms 

dtR^XX) - ^\D X \R 2 (X,Y)) and (d t R 2 (X,Y) - -^-^\Dx\Ri(X, Y)) 

are the second order derivative of analytic functions with a small e prefactor. 
Therefore, one observes that, up to a perturbation, the KH problem behaves 
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like a second order constant coefficient elliptic equation. This fact has several 
important consequence. 

1. It explains why the evolution equation is well-posed only for a short 
time, and with initial data that belongs to the class of analytic functions. It 
is like solving an elliptic equation simultaneously with both the Dirichlet and 
Neumann boundary conditions. 

2. It is a tool for the construction of the solutions that blow up in finite 
time. 

3. It explains, by an indirect regularity argument, the very singular behavior 
of the solution after the first break down of its regularity. 

These three points are discussed in further details below. 

7.2.1 Local solution 

When the curve £t is a graph of a function, say y = y(x, t), the equations (|109| 
and (|110|) become 

Vt + VxVi=v 2l d t uj + d x (viu>) = (123) 
n 1 f y(x,t) -y(x',t)w(x',t) , 

. , 1 f (x — x')u>(x' , t) , , ,„„ . 

where (t>i,i>2) = v is the average velocity given in (|108p . Therefore, the above 
evolution equations involve two unknowns y(x,t) and u>(x, i), which is also the 
case for the Birkhoff-Rott equation (|113[) , where the two unknowns are the two 
components of r(s, t) = (x(s, t),y(s, t)), or of z(s, t) — x(s, t) + iy(s, t). In fact, 
since the Birkhoff-Rott equation has been obtained by choosing the density of 
vorticity as a parameter, one recovers this vorticity by the formula 

, \ 1 



\d s z(s,t)\ 



Since the system is a local perturbation of a second order elliptic equation then 
imposing two constraints at t — is similar to solving this elliptic equation 
simultaneously with both Ncuman and Dirichlet boundary conditions. It is 
known that in the absence of stringent compatibility conditions (they are related 
by the so called Dirichlet to Neumann operator) such a problem can be solved 
only locally and with analytic data. This is the reason why the solution of 
(|123[) - (|125[) is obtained locally in time under the assumption that the functions 
y(x,0) and u>(x,0) are analytic. 



7.2.2 Singularities 

For the construction of singularities one follows the same idea, and furthermore, 
uses the time reversibility of the 2d Euler equations. More precisely, if one 
constructs solutions which are singular at t = and are regular on the interval 
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(0, T] this will imply, just by changing the time variable t into T—t, the existence 
of smooth solutions at t = that blow up at t = T. The first result in this 
direction was obtained by Duchon and Robert [22] ; the initial condition on the 
vorticity at t = is relaxed, and one assumes that the solution y(x,t) goes to 
zero as t — > oo. Then one can consider the system (|123p - (|125[) as a two point 
Dirichlet boundary- value problem with y(x,t) is given for t = and is required 
to tend to zero as t — > oo. Then by a perturbation method one proves the 
following proposition. 

Proposition 7.2 There exists e > such that for any initial data that satisfies 
the estimate 

y(x, 0) = J e^giOdt , with J | 5 (0|d£ < e , (126) 

the problem 1123]) — \125\) has a unique solution which goes to zero as t — > oo . 
Furthermore, this solution is analytic (with respect to (x,t)) for all t > . 

As mentioned above, this is a result about singularity formations. It exhibits 
(by changing the time variable t into T — t) an example of solutions which are 
analytic at some time, but with no more regularity at a later time than what 
is allowed by the equation (|126[) . In fact, it was observed in some numerical 
experiments [53] and [56] that the first break down of regularity appears as a 
cusp on the curve r(A, t) . This motivated Caflsich and Orellana [9] to introduce 
the function 

/o(A, t) = (1 - i){(l - e-i- iX ) 1+ ° - (1 - e-^ x ) 1+ °} , (127) 

which enjoys the following properties 

i) For any t > the mapping A i— ► /(A, t) is analytic. 

ii) For t = 0, the mapping A i— > /(0, A) does not belong to the Holder space 
C 1+(T , but it belongs to every Holder space C 1+a with < a' < a . 

iii) The function 

z {X,t) =X + ef (X,t) 

is an exact solution of the linearized Birkhoff-Rott equation. More precisely, 
one has 

J_ [ f(X,t)-f(X ',t) 
2tt P - V 7 (A - A') 2 



Therefore, by using the cllipticity of this linear operator, one can prove by a 
perturbation method the following proposition. 

Proposition 7.3 For e > 0, small enough, there exist a function r £ (A, t) with 
the following properties: 

i) The function A i— ► r e (X, t) is analytic for t > . 

ii) The function X + e(/o(A, t) + r e (X,t)) is a solution of the Birkhoff-Rott 
equation ftll'ty . 

ii) The function X i— * r e (X,t) is (for X G R, t G R + ) uniformly bounded in 

C 2 . 
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As a consequence of Proposition 17.31 (and of the reversibility in time) one can 
establish the existence of analytic solutions to the Birkhoff-Rott equation (| 1 13[) , 
say in the interval < t < T, such that at time t = T the map A i— > z(X, t) does 
not belongs to C 1+a at the point A = . 

7.2.3 Analyticity and pathologic behavior after the break down of 
regularity 

The local reduction of the KH to the equation 



requires obviously some hypotheses on the regularity of the function z(X,t) near 
the point (0,0). However, when this reduction is valid it will, thanks to the 
ellipticity, imply that the solution is C°°, and even analytic. Therefore, there 
appears to be a threshold (say T) in the behavior of the solutions of the KH. 
Existence of such regularity threshold is common in the study of free boundary 
problems. This threshold is characterized by the fact that any function having 
a regularity stronger than T is in fact analytic, and that there may exist solu- 
tions with less regularity than T. This has the following, practical, important 
consequence: regularity of the solutions that are smooth for t < T and singular 
after the time t = T cannot be extended for t > T by solutions which are more 
regular than the threshold T. Otherwise, the above theorem would lead to a 
contradiction. This fact explains why after the break down of regularity, the 
solution becomes very singular. 

For instance, it was shown by Lebeau [44] (and Kamotski and Lebeau [37] 
for the local version) that any solution that is near a point belongs to Cf(Cj i + ' T ) 
must be analytic. As a consequence, if a solution constructed (by changing the 
variable t into T — t) according to the method of Caflisch and Orellana could 
be continued after time t = T, it would not be in any Holder space C 1+a . 

Therefore, the challenge (and an open problem) is the determination of this 
threshold of regularity that will imply analyticity. Up to now, the best (to the 
best of our knowledge) known result is due to S. Wu [7T] [75]. The hypothesis 
C ; " c (R t ; C]^P{R\)) is replaced by Hl oc (R, t x R-a) • The estimates are done by 
explicitly using theorems of G. David [20] saying that for all chord arc curves 
r '. s i — ► £(s), parameterized by their arc length, the Cauchy integral operator 



is bounded in L 2 (ds) . 

The importance of this improvement is justified by the numerical experiment 
of 39J. It is interesting to notice that these results will apply to logarithmic 
spirals r = e e , 9 6 R, but not to infinite length algebraic spirals. What is 
observed is that from cups singularity the solution evolves into a spiral which 
behaves like an algebraic spiral, and therefore has an infinite length. The results 




dX' 



+E{z(X,t)) (129) 



z(\',t)eu (A - A')(l - e 



/(A,t)-/(A',t) 



A-A' 
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of |7fl| / provide an explanation of the fact that the spiral has to be of an infinite 
length. 

After the appearance of the first singularity the solution becomes very irreg- 
ular. This leads to the issue of the definition of weak solutions (solutions which 
are less regular than the threshold T ) not of the Euler equations themselves, 
but of the Birkhoff-Rott equation. For instance, S. Wu [7T] and [75] proposed 
the following definition: 

A weak solution is a function from R into C,ai-> z(a, t), for which the following 
relation holds 

- Z {a M a) d a) = ± Jj jgg^ *^, 
for every n £ C^°(R) . 

However, the problem is basically open because we have no theorem concern- 
ing the existence of such a solution. Furthermore, for physical reasons weak 
solutions of the Birkhoff-Rott equation should provide weak solutions of the 
incompressible Euler equations, and in fact this is not always the case as it is 
illustrated by (cf. [48]) the Prandtl-Munk example: Start from the vortex sheet 

uo{x 1 ,x 2 ) = —=={x(-i,i)(xi) <8> S(x 2 )) (130) 
v/1 -x{ 

where X(-i,i) is the characteristic function of the interval (—1, 1). By virtue of 
the the Biot-Savart law, the velocity v is constant 

« = (0,~). (131) 
The solution of the Birkhoff-Rott equation is given by the formula 

xx(t)=x 1 (0), x 2 (t) = ~, u)(xi,x 2 ,t)=u (xi,X2 + -). 

On the other hand, it was observed in [35] that the velocity u associated with 
this vorticity is not even a weak solution of the Euler equations. In fact, one 
has 

V ■ u = and d t u + V x ■ (u ® u) + Vp = F , (132) 
where F is given by the formula 

F = £ (Wi + 1, X2 + \) - 6( X1 -1,X2 + t)),0^) . (133) 

This has led Lopes, Nussenweig and Sochet [JS] to propose a weaker definition, 
which contains more freedom with respect to the parameter, and may be more 
adapted. 
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